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1. Introduction 

The experimental realizations |I[ of Bose-Einstein condensation (BEC) |2j in trapped 
dilute alkali atomic vapours have attracted much attention recently. Standard textbooks || 
on statistical mechanics typically treats the topics of non-interacting free Bose gas in one, 
two, and three dimensional boxes in the thermodynamic limit. BEC is defined as the onset 
of a macroscopic occupation of the ground state at a finite temperature, and it is usually 
accompanied by a cusp in the temperature dependence of the specific heat HHJ]. 

Recent theoretical studies on BEC have revealed interesting effects of finite number of 
particles and anisotropic traps fl5 Hl8|l . van Druten and Ketterle @ introduced the concept 



of two-step BEC to explain the observation that the cusp in the specific heat appears at a 
temperature T cusp significantly higher than the BEC critical temperatue T c in anisotropic 
harmonic traps. This phenomena comes about from the changes in the occupations of 
states in the anisotropic trap with different characters as the temperature is lowered. In 
the present work, we study the intermediate steps towards BEC in details. We introduce 
the concept of dividing the set of all single particle states into eight subsets, each of which 
consisting of states with different characters. Such a division is especially useful in treating 
particles in anisotropic traps in that the states in the subsets can be classified into those 
corresponding to zero-point motion, excited states of oscillators in one, two, and three 
dimensions. The occupations in these subsets can be evaluated numerically or by invoking 
a continuum approximation using the concept of the density of states. The possible multi- 
step phenomena associated with BEC can then be described in terms of the changes in the 
relative populations in the different subsets as the temperature is lowered. We found that 
pheneomena involving more than two steps can occur in a general class of highly anisotropic 
traps. We also found that two-step BEC is possible in anisotropic traps with oj\ = uj 2 <C ou 3 , 
in contrast to previous claims ||]. The position of the cusp in the specific heat can also 
be studied within the context of subsets of states. The occupations in the subsets, within 
the continuum approximation, can be expressed in terms of Bose functions g n (z). The Bose 
functions can be classified into two classes according to their behaviour at z — 1. Our 



results indicate that the cusp in the specific heat appears in the vicinity of the crossovers 
in the populations between subsets of states involving Bose functions in different classes as 
temperature decreases. 

The plan of the paper is as follows. The eight subsets of states are introduced in Sec. 2 
together with the method of numerical calculations. Section 3 gives a short discussion 
on the density of states within the continuum approximation. In Sec. 4, we present our 
results demonstrating the possibility of multi-step phenomena towards BEC in anisotropic 
harmonic traps. Section 5 gives a discussion on the cusp in the specific heat in terms of 
the crossover between populations in subsets involving Bose functions of different classes. A 
brief summary is given in Sec. 6. 
2. Eight subsets of single-particle states 

We consider non-interacting bosons trapped in a three-dimensional (3D) anisotropic har- 
monic traps described by the characteristic frequencies uji, cu 2 , and u 3 in the x, y, and 
£-axis. The single particle energy spectrum can be characterized by the quantum num- 
bers rii, n 2 , and n 3 . To facilitate our discussion, we divide the set S = {(ni,n 2 ,n 3 ) : 
ni,n 2 ,n 3 = 0, 1,2, ■••} of all single particle states into eight subsets with different prop- 
erties, i.e., S = S ®S 1 ®S 2 ®S 3 ®S 12 ®S 13 ®S 23 ®S 123 . Here S = {(0,0,0)} con- 
sists only of the overall ground state. The subsets Si = {(ni,0,0) : n\ = 1,2, •■■} , 
S 2 = {(0,n 2 ,0) : n 2 = 1,2, •••}, and S3 = {(0,0, n 3 ) : n 3 = 1,2,---} consist of states in 
which the lowest energy state of the corresponding harmonic oscillators in two of the three 
spatial directions are occupied and the excited states in only one spatial dimension are oc- 
cupied. Hence, Si, S 2 , and S3 include states which are one dimensional (ID) in character. 
Similarly, the subsets S12, S13, and S 23 consist of states which are two dimensional (2D) 
in character in that they consist of states in which one and only one of the three quantum 
numbers vanishes, and hence include states which are excited in the xy, xz, and yz planes, 
respectively. The subset Si 2 3 = {(ni,n 2 ,n 3 ) : ni,n 2 ,n 3 = 1,2,---} consists of all the 3D 
excited states with nonzero rii, n 2 , and n 3 . 

The mean occupation number n p in a single-particle state with energy e p is given, within 
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the grand-canonical ensemble, by 



U P e 0(c P -ri - 1' ^ 



where (3 = 1/ksT is the inverse temperature with ks being the Boltzmann constant. The 
chemical potential \x or the fugacity z = exp(/3/x) is determined by the mean number of 



particles N in the system via iV = J2 P n p- Defining |17j the factor Bj = J2 P exp(— jf3e p ) 



where the summation is over all the single-particle states labelled by p, then we have 

oo 

N=^B d . (2) 

Note that all the effects of the confining potential on the thermodynamic properties are 
included in the factor Bj through the energy spectrum e p . For the 3D parabolic potential 
with possibly different frequencies in different directions, the single-particle eigen-energies 
are given by e ni = Sf=i nihu>i, where uj\ < uj 2 < u 3 are the frequencies characterizing the 
confinement in different directions. The quantum numbers ni, n 2 , and n 3 take on non- 
negative integers, and the zero-point energy has been absorbed into the definition of the 
zero of energy. Once the fugacity z is determined by Eq.(^), the internal energy U can be 
calculated directly via U = — ^ (dBj /d/3) / j, and the heat capacity C v = (dU/dT) N v 
can also be obtained. The subscript v in C v indicates that this quantity is analogous to 
the heat capacity at fixed volume in the sense that the oscillators' frequencies are held 
fixed. The occupation in each of the subsets of states can also be evaluated as a function 
of temperature. Equation (2) forms the basis of our discussion of multi-step phenomena 
associated with BEC. 

3. Continuum approximation using DOS 

In addition to studying the mean occupation in each subset of states numerically, we 
can also study the occupations analytically by invoking the notion of the density of states 
(DOS) P,^,|6HlT| . The DOS for each of the subsets can readily be written down. However, 
one should keep in mind the fact that the different subsets of states have different lower 
limits of energy when a sum over states is turned into an integral, a point that has become 



a source of confusion in the literature 0,0. The mean occupation numbers in the eight 
subsets of states can then be expressed as: 

N w g (z), 



N t « ± gi {ze- a !), £ = 1,2,3, 



(3) 



^« ~ nfefe( ze 0, * = 12, 13, 23, 



Here fij = huJi/ksT (i = 1,2,3), fi* = e*/ksT, and the proper lower limits of energy 
for integrations involving the DOS for subsets corresponding to excited states with ID, 
2D, and 3D characters should be taken to be e* = hui/2 (i = 1,2,3), e* £ = h(ui + u)g)/2 



{it = 12, 13, 23), and e* 23 = h(ui + cj 2 + cu 3 + \Juj\ + u>2+ u;|)/2, respectively. The Bose 
functions g n (z) are defined by g n (z) = Y^jLi z^/j n . These functions can be divided into 
two classes. The first class consists of Bose functions with n > 1, and they are finite at 
z — 1. The second class consists of Bose functions with n < 1, and they diverge at z — 1. 
For example, go(z) = z/(l — z) and gi(z) = — ln(l — z) diverge at z — 1. The results 
in Eq.(3) also give us insight into the multi-step phenomena associated with BEC as the 
dominant subset of states carrying a substantial fraction of the total particles changes from 
one subset to another as the temperature is lowered. Similar approximate expressions for 
the DOS can be obtained in various different ways such as the Euler-Maclaurin summation 
||, the zeta function [[101,1191, or the simple series 0. In the approximation scheme of Ref. 



0, the single-particle states are divided into four groups. 
4. Multi-step BEC 

BEC is usually described as the onset of a macroscopic occupation of the overall ground 
state, i.e., the subset So, at the critical temperature T c . Very often, it is accompanied by the 
presence of a cusp in the specific heat at the same temperature. In recent experiments 
the direct evidence for BEC is taken to be the abrupt increase in the peak density at the 
center of the trapped atomic cloud. Novel experimental techniques have made possible the 
measurements in the spatial and velocity distributions of the trapped bosons. If the subset 
of states containing the largest fraction of particles changes from one subset to another 



consisting of states with a lower dimensional character as temperature decreases, the changes 
in the spatial and velocity distributions should also be measurable. In fact, these changes 
in the characteristic of the trapped bosons have been taken to be the indication of two-step 
phenomena in BEC ||[/]. Here we will study the possible multi-step phenomena in BEC 
due to the changes in the occupations in the subsets of states as temperature decreases Ilj3 



Following Ref.[6], a "saturation" temperature can be defined for each subset of ex- 
cited states. Consider a hypothetical system consisting only of one subset of states, say 
5123, plus the overall ground state. At the saturation temperature T^, the popula- 
tion in the set of excited states saturates. From Eq.(D), is approximately given by 
T 3 £) = (0.94h/k B )(NuJiUJ2 UJ 3) 1 ^ 3 ■ Similarly, the "saturation" temperature for the subsets SV,- 
(ij = 12,13,23) with 2D character is given by T 2 ^> = (0.78h/kB)(Nuj i ujj) 1 ^ 2 ; and that for 
the subsets Si (i = 1,2,3) of ID character is given by T XD ] « (NhoJi/kB)ln(2N) @J^]. For 
example, only the subsets S123 and So are responsible for the BEC in 3D isotropic harmonic 
trap, and the "saturate temperature" T 3 ^ = 0.9AN 1 ^ 3 huj/kB is the corresponding critical 
temperature below which the overall ground state becomes macroscopically occupied. 

For suitably chosen uj x , u> 2 , and u 3 , the competitions among the different subsets may lead 
to complicated and interesting phenomena. Consider 003 : u 2 '■ 0J\ = 2.1 x 10 4 : 4.4 x 10 3 : 1. 
This choice of parameters corresponds to T 3 £> = 1.21T 2D ] and T 2C ] = 2T^. Figure 1 shows 
the typical behavior of three-step phenomena associated with BEC for N = 10 5 particles 
in such anisotropic traps. The temperature is expressed in reduced temperature defined by 
r = kBT/(h 3 cu 1 uj2Uj 3 ) 1 / 3 . As the temperature is lowered, the subset containing the largest 
fraction of particles changes from S'123 to S X2 at about t x = 1.5 (the first step), then from 
S12 to Si at r 2 = 0.6 (the second step), and finally to So at about r 3 = 0.2 (the third step). 
On the other hand, the BEC critical temperature at which N becomes appreciable is at 
about r c = 0.5. Figure 1 also shows the behavior of the heat capacity. The cusp in the heat 
capacity appears at T cusp = 0.65, which is higher than r c . 

Figure 2 shows the results for N = 10 5 bosons in an anisotropic trap with uj x = uj 2 = 
a>3/500 . It was reported in Ref. || that two-step phenomena will not be observed in such 



system since T cusp and r c are nearly the same {r cusp « r c = 0.65). The choice of parameters 
corresponds to the case in which oj% = 002 ^3, and we have T^/T^ = 1.21 (cjg/ 'ufN) 1 ^ m 
1.41. However, our results demonstrated that the particles start to condense from the subset 
S123 into 5*i2 at a temperature t\ ~ 2.0, and then to So at a lower temperature. Such two-step 
phenomena should be observable. 

5. Cusp in specific heat and crossover between subsets 

Our approach of focusing on the occupation in the subsets of states also shed light on the 
position of the cusp in the heat capacity. In the standard textbook case of ideal Bose gas in 
the thermodynamical limit, the temperature T cusp corresponding to the position of the cusp 
coincides with the critical temperature T c of BEC. However, there exists in general no such 
restrictions requiring the coincidence of the two temperatures for systems with finite number 
of particles in anisotropic traps 0. Recall that the Bose functions g n (z) can be classified 
into two classes depending on their behavior at z = 1. For harmonic traps, the occupations 
in 5*123, <Si2, 5*13, 5*23 involve g n (z) with n > 1 (the first class) while the occupations in the 
subsets So, Si, S 2 , S 3 involve g n (z) with n < 1 (the second class). An interesting question is 
that whether there will exist two or more cusps in the specific heat in multi-step phenomena 
associated with BEC. We have studied the temperature dependence of the specific heat for 
a wide range of ratios of ui, u>2, and U3. Our results indicate that there is at most one peak 
in the temperature dependence of the heat capacity, with the peak being rounded off due to 
the finite number of particles in our calculations. The position T cusp of the cusp appears to 
be the results of a combination of crossovers in the occupations between subsets involving 
g n (z) belonging to the two different classes. The crossovers at which the particles condense 
from one subset to another involving Bose functions belonging to the same class do not lead 
to a cusp in the specific heat. For the system considered in Fig.l, one such crossover occurs 
between the subsets S123 and Si at r ~ 0.85, which coincides with a rapid increase in the 
heat capacity as temperature decreases. Another crossover between the subsets S12 and Si 
at r = 0.65 coincides with the cusp in the heat capacity. For the case in Fig. 2, the crossover 
between the two classes of Bose functions occurs between the subsets S123 and So at r ~ 0.65, 
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which coincides with the position of the cusp in the heat capacity. We further re-examine 
the case of 10 6 bosons in a highly anisotropic trap with io 2 = ui^ = 5.6 x lO 6 ^ corresponding 
to Tg^ = 2T2D ■ This system was studied in Ref . || to demonstrate the two-step phenomena. 
Results are shown in Fig. 3. The critical temperature T c is appreciably lower than T cusp . The 
position of the cusp coincides with the range of temperatures in which the two crossovers, 
one between (S12, S13) and Si, and another between S123 and Si, involving different classes 
of Bose functions occur. 

It should be pointed out that our analysis reproduces the standard results in isotropic 
traps. For isotropic harmonic traps in 3D (2D), BEC is accompanied by a crossover in 
occupations between S123 (£12) and S and hence T cusp = T c as S 12 3 (£12) an d So involve 
different classes of g n (z). For ID harmonic traps, BEC is accompanied by the crossover 
between So and Si, both of which involve Bose functions of the same class and hence there 
exists no cusp in the heat capacity. 

Extending our argument to bosons in 3D boxes with dimensions L 1; L 2 and L 3 in different 
directions is interesting . The occupation in a subset of states with (^-dimensional character 
involves the Bose function g d / 2 . Hence only the subset S123 including excited states in all 
three Cartesian coordinates involves g n (z) with n = 3/2 > 1. All the other subsets consist 
of states of lower- dimensional characters and involve g n (z) with n < 1. Thus, BEC occurs 
as long as the spatial dimensions of the box and the number of particles allow the subset 
S123 to be occupied. For a box with Li <C L 2 <C L 3 , we expect a three-step BEC to occur 
with the largest occupation changes, as temperature decreases, from Si 23 to S 23 , then to S 3 , 
and finally to So with the cusp of the specific heat appearing in the vicinity of the crossover 
between S123 and S23. The crossover from S23 to S3 and from S3 to So do not lead to a cusp 
in the specific heat since the subsets involve Bose functions belonging to the same class. For 
a box with L\ <C L 2 = L 3 , we expect a two-step phenomena with the cusp of the specific 
heat appearing at a temperatue near the crossover between the subsets S123 and S23 which 
is higher than the BEC critical temperature. For a box with L\ = L 2 <C L 3 , we expect that 
the cusp to appear when near the crossover between S123 and S3 at a temperature higher 
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than the BEC temperature. 
6. Summary 

In summary, we have studied the intermediate steps towards BEC for finite number of 
bosons confined by anisotropic harmonic traps. Using the concept of dividing the set of all 
single particle states into subsets consisting of states with different characters, two-step and 
multi-step phenomena associated with BEC can be descibed in terms of the crossovers in 
the populations of these subsets as temperature is lowered. Cases with three-step processes 
are demonstrated. The cusp of the specific heat is found to be associated with the crossover 
between subsets involving Bose functions of different classes. The general method adopted 
in the present work can readily be extended to study bosons in other types of traps. 
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FIGURES 

FIG. 1. Numerical results for N = 10 5 particles in a 3D anisotropic harmonic trap with 
ujz : u>2 '■ U\ = 2.1 x 10 4 : 4.4 x 10 3 : 1. The solid lines show the temperature dependence of the 
populations Nq/N in the subset So, Ni/N, N12/N, and N123/N in the subsets S±, S12 and S123, 
respectively. The populations in the other subsets are too small to be shown. The line with circles 
gives the rescaled specific heat as a function of the reduced temperature r = ksT / '(^wiu^^) 1 / 3 . 

FIG. 2. Two-step BEC occurs in a Bose gas of N = 10 5 particles in an anisotropic trap with 
003 : u>2 '■ w\ = 500 : 1 : 1. The particles condense from the subset S123 into the subset £12 at a 
temperature significantly higher than the BEC critical temperature, showing features of a two-step 
BEC. However, the temperature for specifc heat cusp is close to the BEC critical temperature. 

FIG. 3. Typical behaviour of two-step BEC of N = 10 6 particles in a highly anisotropic har- 
monic trap with U3 = L02 — 5.6 x 10 6 oji, corresponding to = 2T^. The condensations from 
the subsets S'123, S12, and Sis to S\ are responsible for the cusp in the specifc heat. 
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